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, Exact solutions to the d-dimensional Schrodinger equation, d > 2, for Coulomb plus harmonic 

oscillator potentials V{r) — —a/r + br^, b > and a 7^ are obtained. The potential V{r) is 
considered both in all space, and under the condition of spherical confinement inside an impenetrable 
^ , spherical box of radius R. With the aid of the asymptotic iteration method, the exact analytic 

■^L ' solutions under certain constraints, and general approximate solutions, are obtained. These exhibit 

the parametric dependence of the eigenenergies on a, b, and R. The wave functions have the 
simple form of a product of a power function, an exponential function, and a polynomial. In order 
to achieve our results the question of determining the polynomial solutions of the second-order 
differential equation 
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Tjj- ■ I. INTRODUCTION 

' A. Formulation of the problem in d dimensions 

The d-dimensional Schrodinger equation, in atomic units fi = /i = 1, with a spherically symmetric potential V{r) can 
be written as 



-l^d + Vir) 



V(r) = E^{r), (1) 



where A^j is the d-dimensional Laplacian operator and — ^1- Following jjyj, in order to transform ([T} to the 

d-dimensional spherical coordinates {r,9\,92, ■ ■ ■ , Od-i), we separate variables using 

V;(r) = T^''-'^'^u{r)Yi,_,..,, {9,... 9d-i), (2) 

where (^i ■ ■ • fi'd-i) is a normalized spherical harmonic with characteristic value l{l + d— 2), / = 0, 1, 2, . . . (the 

angular quantum numbers), one obtains the radial Schrodinger equation as 



1 / (fc-l)(fc-3) 



2 V dr^ 4r2 



V{r) - E 



u(r) = 0, / u'^{r)dr ^l,u{0)^0, (3) 
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where k = d+2l. Assume that the potential V{r) is less singular than the centrifugal term so that 

u(r)^r5(fc— 1) (r — > 0). 
We note that the Hamiltonian and boundary conditions of Q are invariant under the transformation 

{d,l) ^ (d=F2,/±l). 

Thus, given any solution for fixed d and /, we can immediately generate others for different values of d and /. 
Further, the energy is unchanged if k = 2£ + d and the number of nodes n is constant. Repeated application of this 
transformation produces a large collection of states, the only apparent limitation being a lack of interest in some 
values of d (see, for example l2l|). In the present work, we consider the Coulomb plus a harmonic oscillator potential 

V(r) = +br'^, b>0 (4) 
r 

where r = ||r|| denotes the hyper-radius, and the coefficients a and b are both constant. 



B. Degeneracy in spherically confined d-dimensional quantum model systems 

Since the early days of quantum mechanics there has been interest in studying the Schrodinger equation with 
model systems in higher spatial dimensions The so-called accidental degeneracy of the hydrogen atom and 

isotropic harmonic oscillator, characterized by different sets of parity conditions, is generally understood in terms of the 
corresponding 50(4) and SU{3) symmetry groups Following the introduction of 'interdimensional degeneracies' 

flo| there have been several reports involving arbitrary d-dimensional analyses covering many branches of chemical 
physics which have been briefly reviewed in Refs. [ll| - [T^ . It is interesting to note here that the information-theoretical 
uncertainty-like relationships in terms of the Shannon entropy [l5l . [l^ and the Fisher measure [itI [18| are also stated 
in d-dimensional form. 

Owing to the recent interest in quantum dots and fullerine encapsulated electronic systems there has been an 
upsurge of interest in studying model quantum systems confined inside an impenetrable sphere of radius R. We shall 
present here a brief description of the new degeneracy-related changes which are known to occur in the d-dimensional 
H atom Vc = —a/r and the isotropic harmonic oscillator Vh — br^. The eigenspectrum of the spherically confined H 
atom (SCHA) is characterized by three kinds of degeneracy [T^. Two of them are generated from the specific choice of 
the radius of confinement R, chosen exactly at the radial nodes corresponding to the free hydrogen atom (FHA) wave 
functions. In the incidental degeneracy case, the confined (i^, i) state with the principal quantum number ly is iso- 
energic with {i' + 1,£) state of the FHA with energy — l/{2(i/-|- 1)^} atomic units (a.u.), at an R defined by the radial 
node in the FHA. For example, the {v, £) state corresponding to the lowest energy value, when confined at the radius 
R given by the radial node the first excited free state {v +!,£), increases in such a way that the confined-state energy 
becomes the same as excited free-state energy. The specific node in question is given by i? = 0. 24(2^-1- (i—I)(2f-|-(i-|-I). 
Such a degeneracy can be realized at similar choices for R where multiple nodes exist in the second and higher excited 
states of a given £. However, such closed analytical expressions for the radial nodes are not available in the case of 
higher excited states. In the simultaneous-degeneracy case, on the other hand, for all > £ -f 2, each pair of confined 
states denoted by (i^ji) and {n + 1,£ + 2) state, confined at the common R = 0.24(2£ + d — 1){2£ -t- c? -f I), become 
degenerate. Note that the pair of levels in the free state are nondegenerate. Both these degeneracies have been 
shown to result from the Gauss relationship applied at a unique Rc by the confluent hypergeometric functions 
that describe the general solutions of the SCHA problem. Finally, the interdimensional degeneracy [§, [l3| arises, 
as in the case of the free hydrogen atom, due to the invariance of the Schrodinger equation to the transformation 
(£, d) — )■ {£ zt 1, d =F 2). In order to preserve the number of nodes in the radial function, it is simultaneously 
necessary to make the transformation ly — > v -\- 1. The incidental degeneracy observed in the case of a spherically 
confined isotropic harmonic oscillator (SCIHO) is qualitatively similar to that of the SCHA. For example, the only 
radial node in the first excited free state of any given I for d-dimensional SCIHO is located at i? = [2£ + d) /2. For 
the multiple node states, the corresponding numerical values must be used. However, the behavior of the two confined 
states at a common radius of confinement is found to be interestingly different [l^lllj. In particular, for the SCIHO 
the pairs of the confined states defined by {v = £ + 1,£) and [v = £-\-2,£ + 2) at the common R = ^ {2£ + d)/2 a.u., 
display for all v, a constant energy separation of exactly 2 harmonic-oscillator units, 2huj , with the state of higher £ 
corresponding to the lower energy. It is interesting to note that the two confined states at the common R with A£ = 2, 
considered above contain different numbers of radial nodes. The condition for interdimensional degeneracvjol. [l0| due 
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to the invariance of the Schrodinger equation remains the same as before. Recently, the confined systems of the 
d-dimensional hydrogen atom [22| and harmonic oscillator |23| have been studied. Problems involving short-range 
potentials in d dimensions have recently been considered [23.l25{ . In the light of the above discussion, it is interesting 
to study the various aforementioned degeneracies in the free and spherically confined d-dimensional potential generally 
given by V{r) — Vc + Vh = —a/r + br^. 



C. Organization of the paper 



The present paper is organized as follows. In section 2, we discuss some general spectral features and bounds, in 
section 3 we briefiy review the asymptotic iteration method of solving a second-order linear differential equation 
where we discuss the necessary and sufficient conditions for certain classes of differential equations with polynomial 
coefficients to have polynomial solutions. In sections 4 and 5, we use the asymptotic iteration method (AIM) to study 
how the eigenvalues depend on the potential parameters a, 6, R, repectively for the free system {R = oo), and for finite 
R. In each of these sections, the results obtained are of two types: exact analytic results that are valid when certain 
parametric constraints are satisfied, and accurate numerical values for arbitrary sets of potential parameters. 



II. SOME GENERAL SPECTRAL FEATURES AND ANALYTICAL ENERGY BOUNDS 



We shall show shortly that the Hamiltonian H is bounded below. The eigenvalues of H may therefore be char- 
acterized variationally. The eigenvalues Ef^ g = E(a, b, R) are monotonic in each parameter. For a and b, this is a 
direct consequence of the monotonicity of the potential V in these parameters. Indeed, since dV/da = — 1/r < and 
dV/db = r2 > 0, it follows that 

^^%^<0 and ^M^>0. (5) 
oa ob 

The monotonicity with respect to the box size R may be proved by a variational argument. Let us consider two box 
sizes, Ri < i?2 and an angular momentum subspace labelled by a fixed I. We extend the domains of the wave functions 
in the finite-dimensional subspace spanned by the first N radial eigenfunctions for R = Ri so that the new space W 
may be used to study the case R — R2. We do this by defining the extended eigenfunctions so that ipi{r) = for 
Ri < r < i?2- We now look at H in W with box size i?2. The minima of the energy matrix [{ipi, Hipj)] are the exact 
eigenvalues for i?i and, by the Rayleigh-Ritz principle, these values are one-by-one upper bounds to the eigenvalues 
for i?2. Thus, by formal argument we deduce what is perhaps intuitively clear, that the eigenvalues increase as R is 
decreased, that is to say 

OEi^ < 0. (6) 
oR 

From a classical point of view, this Heisenberg-uncertainty effect is perhaps counter intuitive: if we try to squeeze the 
electron into the Coulomb well by reducing i?, the reverse happens; eventually, the eigenvalues become positive and 
arbitrarily large, and less and less affected by the presence of the Coulomb singularity. 

For some of our results we shall consider the system unconstrained by a spherical box, that is to say R = 00. For these 
cases, we shall write E'^^ = E{a,b). If a very special box is now considered, whose size R coincides with any radial 
node of the R — 00 problem, then the two problems share an eigenvalue exactly. This is an example of a very general 
relation which exists between constrained and unconstrained eigensystems, and, indeed, also between two constrained 
systems with different box sizes. 

The generalized Hcisenberg uncertainty relation may be expressed [2^ [l^l for dimension d > 3 as the operator 
inequality —A > (d — 2)^/(4r^). This allows us to construct the following lower energy bound 



E > E ~ min 

0<r<R 



8r^ r 



(7) 



Provided 6 > 0, this lower bound is finite for all a. It also obeys the same scaling and monotonicity laws as E itself 

dr ^ dr 



But the bound is weak. For potentials such as V{r) that satisfy ^(f^^) > 0, Common has shown [28|] for the ground 



state in d = 3 dimensions that (—A) > (l/(2r^)), but the resulting energy lower bound is still weak. 



FIG. 1: The energy E for a — l,b — ^,d — Sasa function of L = ^ for n — 0,1, 2. 



For the unconstrained case R = oo, however, envelope methods j29| - l33l l35j allow one to construct analytical upper 
and lower energy bounds with general forms similar to ([7|). In this case we shall write E^g — E{a, b). Upper and lower 
bounds on the eigenvalues are based on the geometrical fact that V{r) is at once a concave function V{r) = g''^\r^) 
of and a convex function V{r) = g^'^\—l/r) of — 1/r. Thus tangents to the g functions are either shifted scaled 
oscillators above V{r), or shifted scaled atoms below V{r). The resulting energy-bound formulas are given by 



mm 

r>0 



1 

2^ 



a 



b{Pir) 



< Eii{a,b) < min 

r>0 



1 

2r^ 



a 



5(P2r) 



where (Ref. [S^ Eqs.(l.ll) and (1.12a)) 

Pi^n + e+{d-l)/2 and P2 = 2n + e + d/2. 



(8) 



(9) 



We shall sometimes use also the convention of atomic physics in which, even for non-Coulombic central potentials, a 
principal quantum number v is used and defined by 



i^ = n + £+{d- l)/2, 



(10) 



where n = 0, 1, 2, . . . is the number of nodes in the radial wave function. It is clear that the lower energy bound has 
the Coulombic degeneracies, and the upper bound those of the harmonic oscillator. These bounds are very helpful as 
a guide when we seek very accurate numerical estimates for these eigenvalues. 

Another related estimate is given by the 'sum approximation' [33] which is more accurate than ([S]) and is known to 
be a lower energy bound for the bottom Eq^ of each angular-momentum subspace. The estimate is given by 



E^,ia,b)^Si,{a,b)^rmn 

r>0 



(11) 



This energy formula has the attractive spectral interpolation property that it is exact whenever a or 6 is zero. The 
energy bounds ^ and pT|) obey the same scaling and monotonicity laws is those of E'^g{a,b). Because of their 
simplicity they allow one to extract analytical properties of the eigenvalues. For example, in Fig. 1 we show from 
Eq. lfTTj) approximately how the eigenvalue E^g{l, i) depends on £ for n = 0, 1, 2. 

III. THE ASYMPTOTIC ITERATION METHOD AND SOME RELATED RESULTS 

The asymptotic iteration method (AIM) was originally introduced [s^] to investigate the solutions of differential 
equations of the form 



y" ^ Xo{r)y' + so{r)y, (' ^ ^) 



(12) 
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where Ao(r) and so{r) are C°°— differentiable functions. A key feature of this method is to note the invariant structure 
of the right-hand side of (IT^ under further differentiation. Indeed, if we differentiate (|12p with respect to r, we obtain 



y = My + siy 

where Ai = Ag + sq + A§ and si = Sq + sqXq. If we find the second derivative of equation (|12l) . we obtain 

y*^**' = A2y' + S2y 

where A2 = X'l + si + AqAi and S2 ~ s[ + sqXi. Thus, for (n + 1)*'' and {n + 2)"' derivative of n ~ 1,2, 
have 

= Xn-iy' + 



and 



+ SoAn-1. 



respectively, where 

A„ = X'n-i + Sn-l + XoXn-1 and Sn = 

From (fT5)) and ITBl) we have 

A„?;^"^^^ - A„_iy^"+^^ = 6ny where (5„ = A„s„_i - A„_is„ 



(13) 

(14) 
. , we 

(15) 
(16) 
(17) 
(18) 



Clearly, from ([T5| if y, the solution of (|12p . is a polynomial of degree n, then (5„ = 0. Further, if Sn = 0, then Sn' = 
for all n' > n. In an earlier paper [37| we proved the principal theorem of AIM, namely 



Theorem 1 [37]. Given Xq and sq in C°°{a,b), the differential equation has the general solution 



I ^ / f Sn-l(t) 



dt 



C2 + C1 I exp I / (Ao(t) + 2^(r))dr | dt 



if for some n > 



where A„ and Sn are given by jJ?] ). 



Sn — A„s„_i — A„_is„ — 0. 



(19) 



(20) 



Recently, it has been shown [38| that the termination condition ((20|) is necessary and sufficient for the differential 
equation (|12l) to have polynomial- type solutions of degrees at most n, as we may conclude from Eg. dlSp . Thus, using 
Theorem 1, we can now find the necessary and sufficient conditions [30| for the polynomial solutions of the differential 
equation 



(as.or^ + a3,ir^ + a^.^r + 03,3) y" + (02, o?"^ + 02,1^ + ^2,2) y' - {rifir + ti,i) y = 0, 
where a^.j, fc = 3, 2, 1, j — 0, 1, 2, 3 are constants. These conditions are reported in the follow theorem. 



(21) 



Theorem 2 [|39|] Theorem 5]. The second-order linear differential equation 1121]) has a polynomial solution of 
degree n if 



Tifl = n{n - 1) 03,0 + n 02,0, n = 0, 1, 2, . . . , 

provided 03 q + a| q 7^ along with the vanishing of iji + 1) x (n + 1)- determinant A„+i given by 

fio ai rji 
71 Pi "2 ri2 

72 P2 03 m 



(22) 



+1 



7n-2 fin-2 ^n-l rjn-1 
7n-l Pn-l Oln 
7n l^n 
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where all the other entires are zeros and 



= Tis - n[[n - l)a3a + 02,1) 
a„ = ~n{{n - l)a3_2 + 02,2) 
In = Ti,o - {n - - 2)a3,o + 02,0) 

Tin = -n(n + 1)03,3. 



(23) 



Here ti^ is fixed for a given n in the determinant A„-(_i = (the degree of the polynomial solution). The coefficients 
of the polynomial solutions Unif) — X]i=o satisfies the four-term recursive relation 

{i + 2){i + l)a3,3Ci+2 + Hi + l)a3,2 + (i + l)a2,2] Q+i + - l)a3,i + ia2,i - n.i] 

+ [{i - l)(i - 2)03,0 + 1)02,0 - ri,o] c,_i = 0. (24) 



The results of this theorem go beyond the question of finding the polynomial solutions of the second-order linear 
differential equation 



(02,0?'^ + a2.ir + 02,2) y" + (ai,or + ai,i) y' - to,q y = 0. 



(25) 



Indeed Eq.(|25|) has a nontrivial polynomial solution of degree (exactly) ?i e N (the set of nonnegative integers) if, for 
fixed n, 



T-Q.o =="-("•- 1) 02,0 + ai,o, ri = 0,1,2, 



(26) 



provided ai + o^ g 7^ where the polynomials up to a multiplicative constant, may be readily obtained from the 
three-term recurrence relation: 



yn+2 

with the coefficients given by 

((2n + l)a2,o + ai,o)(2(n + 1)02,0 + ai,o) 



yn+i + C„2/„, n > 



(27) 



Bn 

c„ 



("02,0 + 01,0) 

((2n + 1)02,0 + oi,o)(2n(n + 1)02,002,1 + 2(n -t- l)oi,o02,i - 2ai, 102,0 + oi,oOi,i) 



("02,0 + oi,o)(2no2,o + 01,0) 

{n + l)(2(n -I- 1)02,0 + oi,o)((4o2,2o|o - O2,oo|i)n2 -|- (4o2,oOi,o02,2 - oi,oo|i)n + o^ 002,2 ~ oi, 101,002,1 + 02,oOi,i) 



initiated with 



("02,0 + oi,o)(27ia2,o + 01,0) 



2/0 = 1, 2/1 = oi,oa:^ + Oi,i. 



In the next sections, we shall apply the result of theorem 2 to study the possible quasi-exact analytic solutions for the 
rf-dimension Schrodinger equation ^ for unconstrained and constrained Coulomb plus harmonic oscillator potential 
We shall also apply AIM, theorem 1, to obtain accurate approximations for arbitrary potential parameters, again, 
for the unconstrained and constrained c?-dimensional Schrodinger equation 



IV. EXACT AND APPROXIMATE SOLUTIONS FOR UNCONSTRAINED POTENTIAL V{r) 
A. Exact bound-state solutions of a Coulomb plus harmonic oscillator potential in d-dimensions 

In this section, we consider the c?-dimensional Schrodinger equation 



1 fd^ _ (fc- l)(fc-3) 

2 [d^ 



<i{r) = E^iu'!,i{r), 0<r<oo. 



(28) 
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In order to solve this equation by using AIM, the first step is to transform (|28)) into the standard form p^ . To this 
end, we note that the differential equation (j28p has one regular singular point at r = and an irregular singular point 
at r = oo and, since for large r, the harmonic oscillator term dominates, the asymptotic solution of (|28p as r oo is 
Ur^cx, ^ exp(— Y^6/2 r^); meanwhile the indicial equation of pS]) at the regular singular point r = yields 



which is solved by 



s(s-l)-i(fc-l)(fc-3) = 0, 



si = i(3-fc), S2-i(fc-l) 



(29) 



The value of s, in Eq. ((29| . determines the behavior of u^{r) for r — >■ 0, and only s > 1/2 is acceptable, since only in 
this case is the mean value of the kinetic energy finite [40] • Thus, the exact solution of ([28]) may assume the form 



7-2 



l(fe-l) 



exp(- 



/nW, k = d + 3l, 



(30) 



where we note that u'^i{r) 
differential equation for /„(r) as 



j,l{k 1) as r 0. On substituting this ansatz wave function into (|28p . we obtain the 



rf^ir) + {-2r^V2b + k - l) /;(r) + [{2Eii - kV2b) r + 2a] /„(r) 



= 0. 



(31) 



This equation is a special case of the differential equation pip with 03 = 03^1 = 03 3 = 02,1 = 0, 03 2 = 1, 
02,0 = — 2r^\/26, 02,2 = fc — 1, Tifl = —2Ef^i + k\/2b and ri i = —2a. Thus, the necessary condition for the polynomial 
solutions of Ea. ([3T|) is 



Ei, = {2n' + k)^-, 



n' = 0,1,2,.. 



(32) 



and the sufficient condition follows from the vanishing of the tridiagonal determinant A„+i ~ 0, n — 0,1,2, 
namely 

'/3o OLi 
71 /3i Ct2 
72 1^2 as 

7n-2 Pn-2 

7n-l Pn-1 an 
In Pn 

where its entries are expressed in terms of the parameters of Eq. (j3ip by 

/3„ = -2a, Q!„ = -n(n + fc - 2), 7„ = -2(n' - n + l)%/2&, 



(33) 



where n' = n is fixed by the size of the determinant A„+i = and represent the degree of the polynomial solution of 
Eq. pip . We may note that, since the off-diagonal entries Ui and 7^ of the tridiagonal determinant satisfy the identity 

a,;7i > 0, V i = 1,2,..., 

the latent roots of the determinant A„+i are all real and distinct [4l|. Further, we can easily show that the determinant 
(1331) satisfies a three-term recurrence relation 



A, = /3,_iA,_i -7,_ia,^iA,_2, Aq = 1, A_i = 0, ^ = l,2,. 



(34) 



which can be used to compute the determinant A^ (and thus the sufficient conditions), recursively in terms of lower 
order determinants. In this case, however, we must fix n' for each of the sub-determinants used in computing (j34l) . 
For example, in the case of n' = n = 1 (corresponding to a polynomial solution of degree one), we have 
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A, = 



-2a -(fc-l) 
-2\/2fe -2a 



= /3iAi - 7iaiAo=(-2a)(-2a) - (-2\/2fe)(-(fc - l))(l)=4a2 - 2V2b{k - 1), 



that is, the condition of the potential parameters reads 

2a^ - V2b{k - 1) = 0. (35) 
For n' — n — 2 (corresponding to a second-degree polynomial solution) 

-2a -(fc- 1) 
A3 = -4V26 -2a -2k 
-2\/26 -2a 

= (-2a)[(-2a)Ai - (-4V2&)(-(fc - l))Ao]+8/caV^ = 8a(-a2 + 2V26fc - V2b) 

Consequently, we must have 

a{a^ - 2V2bk + V^) = 0. (36) 
In Table I, we give the conditions on the potential parameters to allow for polynomial solutions, from theorem 2. 



—132^2 — 72a2Ai = (— 2a) 



-2a -(fc-1) 
-4V2fe -2a 



- (-2V2fo)(-2fc)(-2a) 



TABLE I: Conditions on the parameters a and 6 for the exact solutions of Ea. H28|) with iJ^j = (2n + k)y^b/2, k = d + 21. 



n 


A„+i 


= 





= 




1 


20^ - 


{k - l)V2b = 


2 


aia^ - 


- {2k - l)V2b) = 


3 


4a4- 


20a^\/26fe - 18&(1 - fe^) = 


4 


a(a - 


- 5V2b{2k + l)a^ + 46(8fc^ -f 8fc - 7)) = 


5 


8a^ - 


140\/26(fc + l)a'^ + 4b(-65 + 518fc + 259P)a^ - 4506\/26(fc - l)(fc + 3)(fc -Hi) = 



It must be clear that although n, the degree of the polynomial solution, it is not necessarily an indication as to the 
number of the zeros of the wave function (node number): further analysis of the roots of fnif) is usually needed to 
compute the zeros of the wavefunction. 

The polynomial solutions fn'{r) — X]"=o '^^^ easily constructed for each n' since, in this case, the coefficients 
Cj satisfy the three-term recurrence relation (see Eg. (124]) ) 



2aQ -f- 2{n' - i + l)V26c,_i . 
c_i=0, Co = l, c,+i = + + , *-0,l,...,n -1, 



(37) 



where n' is the degree of the polynomial solution. When n' = 0, /o(r) = 1. For the first-degree polynomial solution, 
n' = 1. i = 0, we have 



Cl 



2a 



that is 



fi{r) = 1 



2a 



where 2a^ - (k - l)V2h = 0. 



k - 1 

We may further note for a < 0, there is no root of /i(r) = and the un-normalized wave function reads 

2ar 



z,g,(r) =r^(^+2'-i)exp 



2 2 
a r 



d + 2l-l 



1 



d + 2l-l 



a < 



(38) 



(39) 



which represents a ground-state wave function in every subspace labeled by d and I. For a > 0, there is only one root 
of /i and the wave function 
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which represents a first excited-state in each subspace labeled by d and I. The zero of this wave function is located at 

R=^^^ k = d + 2l, a>0. (41) 

2a 



In both cases, a > or a < 0, the exact eigenvalues are given by 

(2d + 4;-l)' efTij""*! 2 



pd - - pd=F2 _ g {d + 2l + 2) ^-^2 _ « 



For second-degree polynomial solution, n' = 2, i = 0, 1, we have for the polynomial solution, ]i(r) = cq + c-\r + C2r^, 
coefficients 

2a , 2(«- - v'2l^(A: - 1)) 

^° = '' ^^ = -fc:rT c, = 

and the polynomial solution then reads 

^-(^^ = ^-(fe3T) + (fe_i)(2fe-i) ' (43) 



where a(a^ — (2fc — l)v26) = from which wc may conclude that — \f2h(k — 1) > 0. Therefore, the wave function 
/2(r) has either two roots or no root based on the value of a > or a < 0, respectively. For a > 0, we have a 
second-excited state wave function 

<(-) = -^'" '--H" 2(2d + 4Z-l) j (^-^T¥3T+ (, + 2/-l)(2.^ + 4/-l) J' «>0' (44) 
which has two zeros at 

2A; - 1 -I- sj2'k - 1 „ 2fc - 1 - \/2A; - 1 , , „, 

i?i = ^-^^ , i?2 = ^-^^ , k = d + 2l, a>0. (45) 

2a 2a 

For a < 0, we have a ground-state wave function 

<^'-> ^A-mr^) ('-^ ^..2,-ix;I^4,-i) )- i^") 

In either case, a > or a < 0, the exact eigenvalues reads 

E^,, ^ Ei ^ Etfi, = -y + ^^ + ^l lini E^,^J, = ^. (47) 
°' ^' 2(2d + 4?-l)' d^oo 2;±i 4 V ; 

For third-degree polynomial solution, n' = 3,i = 0, 1,2, we have for the polynomial coefficients /3(r) = cq + cir + 
C2r^ + c^r^ that 

2a _ (2a^ - 3\/26(fc - 1)) 2a(2a^ - 7^/26fc -|- 3\/2b) 

Co-1, Ci--^— ^, C2- _ 1) ' 3(fc-l)fc(fc + l) ' 

and the polynomial solution then reads 

- 1 - J^r + (2a^-3(fc-l)^) ,2 _ 2a(2a^ - (7fc - 3)^/26) , 

fc-1^ fc(fc-l) 3(fc-l)fc(fc + l) ' ^ ^ 

where the potential parameters satisfy the condition 4a^ — 20a^\/26fc -|- 186(A;^ — 1) = which may be solved in terms 
of \/26 as 

Prom this we have 

f+(\--\ 20,2 (2fc - 3 -h A/T6FT9)r2 4 a^{2Qk'^ - 15fc + 9 + {7k - 3)Vl6fc2 -|- ^)r^ 

■'3 (^-' ~ ~ fc^ " ^ (A: + l)fc(A:-l) ^ 9 3(A; - l)2/fc(A; + 1)2 ' (^°^ 
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and 

2ar (2/fc - 3 - yT6FT9)r2 4 a^{2Qk'^ - IS/c + 9 - (7fc - 3)Vl6fc2 + Q)r^ 

'"^^ ^ (fc + l)fc(fc- 1) ^9 3(fc- l)2fc(fc + 1)2 ■ ^^^^ 

The polynomial f^{r) has two roots if a > and only one root if a < for all r > 0; while f^{r) has no root for 
a < and has three roots for a > (the results that follow from Descartes' rule of signs). In each of these cases, the 
eigenvalues are given by 



2 (d + 2/ + 6) ( 5d + 10/ ± ^16(d + 2Z)2 + 9 



^" = ¥- (d + 2/-l)(. + 2/ + l) ^' 

We can also show for the fourth-degree polynomial solution, rt' = 4, i = 0, 1, 2, 3, we have 

2a (2a2 - 4V26(fc - 1)) 2 4a[-a2 + 5%/2fefc - 2\/26] 3 2[a4 - (1 + 8fc)%/26a2 + 125(fc2 _ 1)] ^ 
•''^^''^ ~ ^ ^ V^'' ^ fc(fc- 1) 3(fc- l)fc(fc+ 1) ^ 3{k + 2){k + l)k{k - 1) 

(53) 

subject to a{a^ - 5%/26(2fc + l)a2 + 46(8fc2 + 8fc - 7)) = and in this case 

^2 (fc + 8) (lOfc + 5 ± 3v/(2fc+ 1)2 +8) 
^"^y 2(2fc + l)2^9 ^' ^-'^ + 2/,a^0 (54) 

and similarly for other cases. Indeed, using the recurrence relation, Eq.(j37l), it is straightforward to compute explicitly 
the polynomial solution of any required degree. 



B. Approximate solutions for arbitrary potential parameters on half-line 

For arbitrary values of the potential parameters a and b that do not necessarily obey the above conditions, we may 
use AIM directly to compute the eigenvalues accurately, as the zeros of the termination condition (^0)) . The method 
can be used, as well, to test the exact solutions we obtained in the above section. To utilize AIM, we start with 

(55) 

.Soir)^~{2Ei,-kV2b)~^ 

and computing the AIM sequences A„ and s„ as given by Ea. ([T7| . We should note that for given values of the 
potential parameters a, b, and of A; = d + 2l, the termination condition (5„ — A„s„„i — A„_is„ = yields an expression 
that depends on both r and E. In order to use AIM as an approximation technique for computing the eigenvalues 
E we need to feed AIM with an initial value of r = ro that could stabilize AIM (that is, to avoid oscillations). For 
our calculations, we have found that ro = 3 stablizes AIM and allows us to compute the eigenvalues for arbitrary 
k = d + 21 and n as shown in Table HIl There is no magical assertion about rp = 3, indeed using an exact solvable 
case, say E = 2.5 with d = 3, / = 0, n' = 1 for a = 1 and b = 1/2, we may approximate r = ro by means of 
E — V{r) — which yields ro ~ 2.4 as an initial starting value for the AIM process. The eigenvalue computations 
in Table |ll] were done using Maple version 13 running on an IBM architecture personal computer, where we used a 
high-precision environment. In order to accelerate our computation we have written our own code for a root- finding 
algorithm instead of using the default procedure Solve of Maple 13. 



V. EXACT AND APPROXIMATE SOLUTIONS FOR CONSTRAINED POTENTIAL 



A. Analytic solutions 



We now consider the d-dimensional Schrodinger equation 

(fc- l)(fc-3) 



4r2 



Vir) 



uii{r)^Eiiuii{r), < r < R, 



(56) 
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TABLE II: Eigenvalues ^^-2, 3,4,5, 6,7 y^^-^ + /2. The initial value utilize AIM is ro = 3. The subscript N refer to 

the number of iteration used by AIM. 



n 


-C'nO 


Tl 


7Titi = 3 

-C'nO 





-1.836 207 439 051 476 488jv 


=78 





0.179 668 484 653 553 873 jv= 


71 


1 


1.576 895 542 024 474 773jv 


=71 


1 


2.500 000 000 000 000 000iv= 


3 Exact 


2 


3.828 388 290 161 145 035jv 


=64 


2 


4.631 952 408 873 053 214jv= 


59 


3 


5.963 137 645 125 787 098jv 


=60 


3 


6.712 595 725 661 429 760jv= 


58 


4 


8.052 626 115 348 259 660jv 


=59 


4 


8.769 519 600 328 899 714jv= 


57 


5 


10.118 396 975 257 306 974jv 


=58 


5 


10.812 924 292 726 383 736iv= 


=56 


6 


12.169 728 962 611 565 630jv 


=58 


6 


12.847 666 480 105 796 414jv= 


=55 


n 


rpd^i — 
-C'nO = ^nl 


n 


T^d—b — 773 
^nO = J^nl 





1.039 629 453 693 666 062]v= 


63 





1.709 018 091 123 552 219jv= 


60 


1 


3.191 127 807 756 594 984]v^ 


58 


1 


3.801 929 609 626 278 046jv= 


55 


2 


5.273 870 062 099 308 315]v= 


57 


2 


5.860 357 172 819 176 603iv= 


55 


3 


7.329 588 502 119 331 779n= 


54 


3 


7.902 317 748 608 790 676 iv= 


52 


4 


9.371 002 235 676 830 145]v= 


53 


4 


9.934 707 216 855 127 521jv= 


51 


5 


11.403 631 794 919 620 038iv= 


=53 


5 


11.960 878 210 587 747 585iv= 


=51 


6 


13.430 355 223 285 870 950iv. 


=53 


6 


13.982 705 338 927 982 296jv= 


=50 


n 


rpd=6 — ipi — f2 
-C/nO = -t^nl = ^n2 


n 


rpd—7 — jTi5 — 7Ti3 
-C'nO = -C'nl = -C'n2 





2.311 633 609 259 797 633]v= 


56 





2.882 228 025 698 769 118jv= 


52 


1 


4.376 287 059 247 773 643]v= 


52 


1 


4.930 673 420 047 524 772jv= 


50 


2 


6.420 575 455 465 976 803]v= 


52 


2 


6.965 837 318 124 071 248iv= 


48 


3 


8.453 864 208 404 376 310n= 


49 


3 


8.993 183 541 000 280 301iv= 


47 


4 


10.480 309 010 248 775 417iv 


=50 


4 


11.015 405 735 685 306 483iv 


=47 


5 


12.502 107 717 572 917 269jv. 


=48 


5 


13.034 024 645 001 717 876jv= 


=47 


6 


14.520 559 118 487 876 168iv= 


=48 


6 


15.049 980 236 362 344 263iv= 


=47 



where 



br^, iiO<r<R 



Vir) 



(57) 



00 if r > _R 

and u^;(0) = uf^iiR) — 0- We may assume the following ansatz for the wave function 



(i? — r) exp 



^rA /„(r), k = d + 2l. 



(58) 



where R is the radius of confinement, and the {R — r) factor ensures that the radial wavefunction uj^;(r) vanishes at 
the boundary r = R. On substituting (|49|) into (|47| . we obtain the following second-order differential equation for 
the functions /„(r): 



k - 1 



1 



2fo r f'^{r) 



2r R-r 

{-2Ef^i + (fc + 2)V2by + {R{2Ef^i - kV2b) - 2a)r -k + l + 2Ra 



r{R - r) 

We note that this equation reduces to Ea. (l5T|) as i? -> 00. Equation ([55]) can be written as 
[-r^ + Rr]fl[{r) + [2V2br^ - 2V2bRr^ - (fc + l)r + (fc - l)R]f^{r) 

{-2Efa + (fc + 2)V2b)r^ + {R{2E'^i - kV2b) - 2a)r + 2Ra-k + l 



Ur). 



(59) 



/40 =0 



(60) 
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This differential equation cannot be studied using Theorem 2. Consequently a further investigation of the following 
class of differential equations 

(04,0'''' + 04, ir^ + 04, 2?^^ + 04,3?^ + 04,4)2;" + (03,0?"^ + 03,1?'^ + 03,2?^ + 03,3)?/' - (t2,o?'^ + T2,ir + T2.2)y = 0, (61) 

is needed. Indeed, by using Theorem 1 and a proof along the lines of the proof of Theorem 2, we are able to establish 
the following: 

Theorem 3. The second- order linear differential equation i61\) has a polynomial solution y{r) = X]fc=o '^fc'''^ */ 

T2,o 'T-l" - 1) 04,0 + 03,0, n = 0,1,2,..., (62) 
provided 04 q + o| q 7^ where the polynomial coefficients Cn satisfy the five-term recurrence relation 

{{n ~ 2){n - 3)04,0 + (n - 2)03,0 - ■^2.o)c,i-2 + {{n - l){n - 2)04,1 + {n - 1)03,1 - T2,i)c„_i 

+ {n{n - 1)04,2 + ^^03, 2 - T2,2)c„ + (ri(n + 1)04,3 + {n + l)o3,3)c„+i + (n + 2){n + l)o4,4C„+2 = (63) 

with c_2 = c_i = 0. 

In particular, for the zero-degree polynomials co 7^ and c„ =0, n > 1, we have 

T2,2 = 0, T2,i = 0, T2,0 = 0. 

For the first-degree polynomial solution, co 7^ 0, ci 7^ and c„ = 0, n > 2, we must have 

T2fl = 03,0 

along with the vanishing of the two 2 x 2 -determinants, simultaneously, 



'T2,2 03,3 
-T2,l 03,2 — T2,2 



0, and 



'T2,2 03,3 

-03,0 03,1 — r2,i 



0. 



For the second-degree polynomial solution, cq 7^ 0, ci 7^ 0, C2 7^ and Cn = for n > 3, we must have 

T2,0 — 2 04,0 + 2o3,o 

along with the vanishing of the two 3 x ^-determinants, simultaneously, 



— T2,2 03,3 204,4 

—T2,i 03,2 — r2, 2 204,3 + 203,3 
-2o4,o — 2o3,o 03,1 - r2,i 2o4,2 + 203,2 — r2,2 



0, and 



-T2,2 03,3 204,4 

'^2,1 03,2 — r2, 2 204,3 + 203,3 

-2o4,o - 03,0 2o4,i + 203,1 - T2,i 



(64) 
(65) 

(66) 

(67) 

0, (68) 



and so on, for higher- order polynomial solutions. The vanishing of these determinants can be regarded as the conditions 
under which the coefficients T2.1 and T2^2 of Eg. i fgT)] are determined. 



Using Theorem 3, we may note, with 04,0 — 04,1 = 04. 4 = 0,04,2 = —1,04,3 = R,a3fi = 2\/26, 03,1 — — 2'\/2&i?, 03,2 — 
-(fc + 1),03,3 = (fc- l)i?,T2,o = 2E^i - (fc + 2)V26, r2,i = -{R{2Ef^i - kV2b) - 2o),T2,2 = -2i?o + fc - 1, that the 
necessary condition for polynomial solutions fn{r) = X]fe=o '^^'''^ Ea. (j60p is 



Ef^i = ^{2n + k + 2)V2b, k = d + 2l, 



(69) 



where n refers to the degree of the polynomial solution and not necessarily to the number of zeros for the exact wave 
function. For sufficient conditions, we have for the zero-degree polynomial solution n — 0, that Eq. (|64p yields 



/o(r) = l, E^i ^ ^{k 2)V2b, if a ^ RV2b and Ra = ^(k - 1) 



(70) 



where, again, k = d + 21. For example, if o = 3, b — 4.5, we have R = I and for fc = d + 2/ = 7, we have the exact 
solution 



^00 — ^01 — ^02 



13.5 



13 



and for a = 4, & = 8, we have R = 1 and for fc = d + 2Z = 9, we have the exact solution 



Thus, for the values of the potential parameters a, b and R as given by 

{a,b, R) 

we have the exact solutions 



^{2l + d-l),^{2l + d~lf,R 



(71) 



E, 



01 iR 



l,{d + 2l-l){d + 2l + 2), 



<(r) =r^(2'+f'-i)(i?-r)exp(- 



(72) 



d+2l-l Jl 



We may note that the confinement size R — (k — l)/(2a) represent the root of the unconfined wave function, (|40p . 
with the same energy (compare (|42] with (ITO)) ). 



(73) 



For first-degree polynomial solution n = 1, we have using (j69p . or r2.o ~ 4:V2b, 

Ei = l{k + 4)V2b, 

along with the two conditions, obtained using (j66p . which relate the potential parameters by 

2kRa - k{k - 1) - 2i?2a2 + 2R^V2b{k - 1) = 0, 

2\/26i?2 - 2Ra + fc - 1 = 0. 
where, in this case, the polynomial solution reads 

{2Ra+l-k) 
/i(r) = 1 — r. 



R{k-1) 



Thus, for the relations 



^ u u^ /2fc-l + V2fc^ (fc + V2fc^)2 
(a, 6, i?) = ( — , ^7^3 ,R 



2R 



k = d + 2l 



we have the exact solutions 



' ^1/ = 4^2 (A: + 4)(fc + V2fc^), 



uf,(r) =r^('=-i)(i?-r)exp(-^±^ 



and for 



{a,b,R) 



2k -1- V2fc - 1 (fc - V2fc- 1)^ 



2i? ' 8i?4 

we have the exact solutions 

f = ^(fc + 4)(fc- V2A^), 



-, _ (fc+V2FT) 



i? , k = d + 2l, 



^^^(r) ^r^(fc-i)(i?~r)exp(- ^-y^ ^2^ _ (W2feHi). 



(74) 



(75) 



(76) 



(77) 



(78) 



(79) 



We note that these exact-solutions cases ((76)) and (|78|) represent the nodes of the wavefunction in the infinite case 

da. 
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For second-degree polynomial solutions n = 2, we have the exact eigenvalues 

Et = \{k + (;>)V2h (80) 

where k = d + 21 and the potential parameters a, b and R are related by the following two conditions (obtained from 
the two determinants in 



m^a^ - 6(fc + l)i?V - 2R{R'^V2b{7k - 3) - 3k{k + l))a + 3{k - l){k + l){3V2bR'^ - k) ^ 0, (81) 

and 

2R^a^ - 2R{R'^V2b + k)a^ - {k - l){3V2bR^ - k)a + 6b{k - 1)R^ = 0. (82) 

In this case the exact solution reads 

d ^ f 2W1 2Ra-k + l (2R^a^ -2Rak + kik-l)-3V2bR^k~l)) A 

^ni?-Oexp(^-^-.^j [l~ r + j^,^^^ r^j. (83) 

Again in this case we can show that these exact solutions correspond to the zeros of the wavefunction in the infinite 
case (EOl and (El 



Similar results can be obtained for higher n (the degree of the polynomial solutions). It is important to note that 
the conditions reported here are for the mixed potential V{r) = —a/r + br"^, where a ^ and 6 7^ (that is to say, 
neither coefficient is zero). 

B. Approximate solutions for confined potential with arbitrary parameters 

For the arbitrary values of a, b and R, not necessarily satisfying the above conditions, we may use AIM directly to 
compute the eigenvalues with a very high degree of accuracy. This also allows us to verify the exact solutions we 
obtained in the perevious sections. Similarly to the unconfined case, we start the iteration of the AIM sequence A„ 
and Sn with 



AoM = -2te-^-^/26: 



/ N _ (-2Ef^, + {k+2)%/2b)r'^+(R(2E^i-kV2b)-2a)r+2Ra-k+l 
S0(r) - riR-r) ■ 



(84) 



where < r < R. It is interesting to note in this case, that, unlike the unconfined case, the roots of the termination 
condition (5„ = A„s„_i — A„_is„ — are much easier to handle in the present case. This is due to the fact that vq 
is now bound within (0,i?) for every given R. Thus, it is sufficient to start our iteration process with initial value 
To = i?/2. In table Hill we reported the eigenvalues we have computed using AIM for a fixed radius of confinement 
R = 1, with To ~ 0.5 as an initial value to seed the AIM process. In general, the computation of the eigenvalues is 
fast, as is illustrated by the small number of iteration N in Tables Hill The same procedure can be applied to compute 
the eigenvalues for other values of a, b i?, and arbitrary dimension d. The results of AIM may be obtained to any 
degree of precision, although we have reported our results for only the first eighteen decimal places. It is clear from 
the table that our results confirm the invariance of the eigenvalues under the transformation (d, Z) — > (d =F 2, Z ± 1). 



VI. CONCLUSION 



We study a model atom-like system — — a/r which is confined softly by the inclusion of a harmonic-oscillator 
potential term b and possibly also by the presence of an impenetrable spherical box of radius R. For 5 > or i? < 00, 
the entire spectrum i?^ g{a,b, R) is discrete. We have studied these eigenvalues and we present an approximate spectral 
formula for the 'free' case, R = 00. For the general case of i? < 00, AIM has been used to provide both a large number 
of exact analytical solutions, valid for certain special choices of the parameters {a, &, i?}, and also very accurate 
numerical eigenvalues for arbitrary parametric data. In the cases where we have found analytic solutions for i? = 00, 
the exact wave functions are no longer expressed in terms of known special functions, as is possible for the hydrogen 
atom. However, the exact solutions we have found for confining potentials correspond to confinement at the zeros 
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TABLE III: Eigenvalues E'^Y^''^{a, b; R) for V{r) = -a/r + br'^,r G (0, 7?), where a = ±1, b = 0.5, i? = 1 and different n and I. 
The subscript N refers to the number of iteration used by AIM. 



n 


/ 


i:;:r(ia/2:i) 


/) 


/ 










-1.275 615 599 206 285 795n-- 


=32 








- 1.275 615 599 206 285 795jv^a2 




1 


5.400 467 192 272 980 536 jv= 


=26 


1 




10.924 630 155 130 440 587jv=27 




2 


11.652 661 600 597 110 050Ar= 


=24 


2 




32.734 045 433 763 800 052;v=32 




3 


19.010 259 174 813 201 428Ar= 


=24 


3 




64.522 506 980 951 712 401jv=4o 




4 


27.565 689 679 299 850 255jv= 


=25 


4 




106.243 804 852 673 032 613jv=47 




5 


37.324 795 658 776 520 956iv= 


=28 


5 




157.875 359 994 443 580 341jv=53 


n 


/ 


i5^r'(-l,l/2;l) 


n 


I 


B^Jr' (-1,1/2;!) 








6.107 045 323 129 696 121jv= 


=27 








6.107 045 323 129 696 121jv=27 




1 


9.530 081 242 027 809 913jv= 


=24 


1 




19.534 700 629 074 427 546jv=27 




2 


15.106 527 319 660 138 719Ar= 


=24 


2 




42.295 175 016 376 479 090;v=35 




3 


22.149 694 772 638 116 456 Ar= 


=24 


3 




74.728 314 736 027 030 722;v=43 




4 


30.518 339 762 183 359 381iv. 


=27 


4 




116.935 489 978 445 435 860jv=47 




5 


40.151 787 835 702 316 786Ar= 


=29 


5 




168.956 853 183 684 793 313jv=55 


n 


/ 


S^r'(l,l/2;1) 


n 


I 


i5;lr'(l,l/2;l) 








5.400 467 192 272 980 536 jv= 


=26 








5.400 467 192 272 980 536Ar=24 




1 


11.652 661 600 597 110 050iv. 


=24 


1 




22.123 225 647 087 677 088iv=25 




2 


19.010 259 174 813 201 428Ar= 


=24 


2 




48.910 542 938 654 909 374jv=35 




3 


27.565 689 679 299 850 255jv= 


=25 


3 




85.660 358 190 161 408 159iv=43 




4 


37.324 795 658 776 520 956 Ar= 


=28 


4 




132.333 925 295 766 891 686jv=48 




5 


48.278 874 241 139 597 779jv= 


=31 


5 




188.912 601 544 108 537 448iv=64 


n 


/ 


i5^5r"(-l,l/2;l) 


n 


I 


EiTH-hl/2;l) 








9.530 081 242 027 809 913jv= 


=24 








9.530 081 242 027 809 913Ar=24 




1 


15.106 527 319 660 138 719Ar= 


=24 


1 




27.374 386 080 371 192 265iv=27 




2 


22.149 694 772 638 116 456iv= 


=24 


2 




54.884 084 396 689 521 442jv=36 




3 


30.518 339 762 183 359 381iv= 


=27 


3 




92.165 200 694 649 737 766jv=44 




4 


40.151 787 835 702 316 786jv= 


=29 


4 




139.258 753 086 612 471 603jv=49 




5 


51.014 646 696 330 218 668jv= 


=32 


5 




196.184 615 317 703 801 052jv=64 



of the unconfined case. An interesting qualitative feature seems to be that E^i{a,b,R), for large R, is concave with 
respect to n, i, or d, but becomes convex as R is reduced; this may arise because the reduction in R perturbs the 
higher states more severely since, when free, they are naturally more spread out. It is hoped that the work reported 
in the present paper will provide a useful addition to the growing body of results concerning the spectra of confined 
atomic systems in d dimensions. 
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